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JACK MORAVA 

Abstract. The classical Einstein-Hilbert action in general relativity 
extends naturally to a blow-up (in the sense of algebraic geometry) of 
the usual space of pseudo-Riemanriian metrics; this presents the metric 
tensor gut as a kind of Goldstone boson associated to the real scalar field 
defined by its determinant. This seems to be quite compatible with the 
Higgs mechanism in the standard model of particle physics. 



0.1 Gauge theories of mass can be traced back to Weyl's early work on con- 
formal geometry [10 §1, 15 §16], in light of de Broglie's quantum-mechanical 
relation 

mc^ = hv 

between energy and frequency; since then, work in particle physics culminat- 
ing in the standard model has made the subject predictive. This note argues 
that relations between physical mass and (Lorentz-Minkowski) geometry are 
best understood by thinking of the usual pseudometric tensor as the prod- 
uct of a 'dimensionless' unimodular tensor (roughly \<\ei g\~^/^gik\) with a 
suitable power of a density 7 ~ | det g\~^^'^'^ (a section of a real line bundle 
which transforms like an inverse length). This enlarges the class of possible 
states by allowing | det g\ — t- or 00, provided the conformal behavior is 
otherwise good. 

Models of this sort have a long history in physics, but related recent devel- 
opments in Riemannian geometry seem not to have had much impact there. 
In particular, Yamabe's nonlinear elliptic eigenvalue equation [11, 16] 



A + i^i?(5)-AH^/("-^) 
n — 1 







whose solutions u > define conformal deformations g := '^^g of the 

metric g with scalar curvature 

Tl — 1 

Rig) = 4 A = constant 

n — 2 

is (apart from the change from Riemannian to Minkowski signature) iden- 
tical with a (super-renormalizable, if n = 3, 4 or 6) conformally invariant 
nonlinear wave equation with its own distinguished literature (cf. eg [3, 9, 
13 §15.2]). 
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0.2 Aside from an example in §4.2, there is little in the present account of 
this material that is new, but the literature of these questions is old, broad, 
and perhaps confusing. My sense is that the subject is simpler than one 
might think; but there are nontrivial conceptual issues at stake, which are 
hard to present clearly. 

I've tried to simplify things by keeping the language, and especially the 
physics, as classical as possible: in particular, though there is some dis- 
cussion of the Higgs mechanism in §2.2, it occurs here only as a classical 
field theory; delicate quantum-mechanical issues involved in its renormal- 
ization are completely ignored. However, there are some mathematical sub- 
tleties: the nonlinear eigenvalue problems considered here may admit weak 
(ie moderately non-smooth) solutions, which correspond to phase transi- 
tions in some physical models. The final example suggests the possibility 
that a change in conformal structure in the interior Schwarzschild region 
might correspond to a change in the physical vacuum state there. 

The paper is organized as follows: §1 summarizes some background, eg a 
blow-up of spaces of quadratic forms at 0, basic facts about densities on 
manifolds, and associated moduli spaces of geometric data. The next sec- 
tion uses this formalism to identify the Einstein-Hilbert action of classical 
relativity with a version of Yamabe's conformally invariant functional. §3 
notes some consequences of the Sobolev embedding C /^^n/(n-2) ^j^j^j^ 
seem to fit with old ideas from physics about Mach's principle, and §4 dis- 
cusses some classical examples in this framework. 

0.3 I'd like to thank S. Agarwala and J. Baez for helpful, patient, and 

skeptical conversations in the early stages of this work, and Dan Christensen 
for unflappable help with confusing numerical matters toward the end. 

§1 Desingularizing quadratic forms at 

1.1 A nondegenerate symmetric bilinear form q of signature /c on a finite- 
dimensional real vector space V has a group 0{q) C Gl(y) of linear isomor- 
phisms. If n is the dimension of V, the homogeneous space Gl(F)/0(g) is 
isomorphic to the space Qni^) of nondegenerate quadratic forms of signa- 
ture k on V . This is an open cone, and its quotient 

Gm\G\{V)/0{q) :- \Q\i{V) 

by the center of Gl(y) is a space of 'unimodular' quadratic forms of signature 
k, isomorphic to the (contractible) space of real n x n symmetric matrices 
whose eigenvalue configuration A* consists of n± = ^(n it k) positive (resp. 
negative) real numbers satisfying 

(_l)(n-fe)/2-Q^^ = 1. 
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For any real s, let be the one-dimensional real representation of Gl^; 
defined by 

sgn(dct) • I dct Y'"" : G1„(M) ^ ; 
it has a complex analog when s G C. 

The quotient Q* of Q* x M(^n — 1) defined by the action 

p-{q, r) ^ (|p|^/("-2)g, p-\) 

of /9 G is a blowup 

(the right-hand map sends [q,r] to (| detg|^^/" • q, \ det g|^"~"^^/^'*r)) which 
replaces the cone of quadratic forms by the corresponding cylinder. 

1.2 For simplicity I'll assume M compact closed, with principal frame bundle 
Gl(rM), and I will write M(s) for the real line bundle 

^{s) Xgi„(m) G1(Tm) ^ M 
with fiber R(s). There is a canonical isomorphism 

A"r^ ^ M(n) 

between n-forms on M and sections of M(n), and hence a canonical Lebesgue 
functional on the space rM.(n) of sections over M; sections of M(0), on the 
other hand, are ordinary real-valued functions. When n > s > 0, let L(s) 
denote the real Banach space of (equivalence classes of) sections of M(s), 
completed with respect to the norm 

ii'^iis := [ / i<^r"/^l'^" 

(or by the essential supremum when s = 0). Holder's inequalities define a 
canonical isomorphism 

L(s)* ^L(n-s) 

and continuous pairings 

L(s) (g)]R L(t) ^ L(s + 1) 

when {) < s,t,s + t < n. rM(s) C L(s) will denote the subspace of smooth 
sections. 

1.3 A pseudometric of signature A; on M is a smooth section of the bundle 

Qi{M) := G1(Tm) Xgi„(r) ^ M ; 
it defines a smooth section 

*^/"(l) := sgn(deg5) • |det5|"/2n g L(g) 

of M(s), as well as an isometry 

(t> ^ = / : L(s) ^ L"/^(M, dvolg) . 
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A pseudometric g and a density (f) G L(|n — 1) define, at each point x G M, 
an element [g{x),(j){x)] S Qn{TM,x) and thus a map 

TQi X rM(in - 1) ^ rg; 

which sends the pair {g, (p) to a section of the bundle of blown- up pseudo- 
metrics. However it is really the ray defined by (p that is significant, and it 
will be useful to assume that 

(p G r'R(in - 1) := m(in - 1) - {0} . 

This constructs a kind of coarse moduli space of generalized pseudometrics, 
which has a finer analog: an element u of the group rK.(0)^ acts on (g, <f)), 
sending it to (|n|^/("~^)(7, n~^(/)), and the map to FQ* factors through the 
quotient of this action. 

§2 A conformally invariant model 

"There is nothing in the world bigger than the tip of an 
autumn hair, and Mount T'ai is little." 

Chuang Tzu, Discussion on making all things equal, tr. 

Burton Watson 



2.1.1 Proposition: The diagram 

rg; x r'M^n - 1) — -rg; x™(o)x t'^'' 

igX7 

FQl 

commutes; where 




2^ 



M 



R{g) dvolg 



is the Einstein-Hilbert action functional (with k = SttG), 



Y[g, 



M 



\df\l + l^R{9)f 



' n 



dvoL 



is Yamabe's conformally invariant quadratic forrr0(with / as in §1.3), and 



n 



Proof: This is an absolute triviality (except for the assertion that Y is 
conformally invariant, which is now classical [11, 16]). Note however that 
the dotted arrow is not asserted to exist. □ 

^Here Yis normalized as if it were the Lagrangian in a Feynman measure of the form 
exp(-iL 

matter 
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2.1.2 This reformulates the Einstein-Hilbert action as the Lagrangian for a 
conformally invariant physical theory involving a unimodular pseudometric 
(ie a section of the bundle ^) and a real spin-zero gauge field 7, defined 
locally by measurements of Newton's constant, with symmetry broken on the 
locus where 7 — > (or 00, if we allow noncompact M). Away from this set, 
conformal invariance allows us to assume that (the function corresponding 
to the density) 7 is constant, ie roughly the Planck frequency 

(|G/i)-^/2 ~ 90.7 X 10^^ MHz . 

This, after all, is what a gauge theory does; we understand 7 to be constant 
because it defines the local mass scale. 

At first sight Y looks like the Lagrangian for a real scalar boson, moving in 
a potential field of the form R{g) [9]; but requiring that (p not be identically 
zero can be interpreted as the introduction of a self-interaction term. The 
Sobolev embedding theorem says that (on a compact smooth n-manifold) 
the space Lf of functions with s derivatives in embeds in iS t — n/q < 

s — n/p: in particular, L\ C Lq"^^"" is just on the edge of continuity. 
Requiring that 7 have fixed norm as a {\n~ l)-density, ie that its Lebesgue 
2n/(n— 2)-norm be finite, is equivalent to adding a Lagrange multiplier term 
of the form 

^(||^||2n/(n-2) _ ^) ^ 

to the Lagrangian; which, when n = 4, is equivalent to allowing the 'dilaton' 
7 a quartic (super-renormalizable) self-interaction. 

2.1.3 From this point of view, the 'graviton' (ie, the field represented by the 
rank two symmetric tensor ^jj.) is a Goldstone boson associated to 7: if we 
'decouple' the metric from its determinant by writing 

with 

(t>={\ det ^|V2)(n-2)/2n ^ p^^l^ _ 

(so I det^l^/^ = 1), then Yamabe's equation 



j *,R{g) = I [<p'R{g) + 4^ 



2 



is completely analogous to Goldstone's identity 

\d{e^\)\l = x'ml + \dx\l . 

The opposite interpretation - that the dilaton is a Goldstone boson as- 
sociated to the metric - is more usual in physics [12]. The interpreta- 
tion here is that the boson associated to the conformally invariant wave 
equation is more significant locally, while (perturbations of) the Lorentz- 
Minkowski metric, though fundamental for geometry, become important 
only at quantum-mechanically vast distances. 
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In fact such issues go back to the earliest days of the subject. Weyl observed 
[15 §28] that the Einstein-Hilbert action can be written as a quadratic func- 
tional 

s{g) := J g^'in^r^, - r^,,r*,,] dvoi, 

in the first derivatives of g, analogous to the left-hand side of Goldstone's 
identity. 

2.2.1 It is well-known, but perhaps quite remarkable, that the standard 
model of particle physics is very close to conformally invariant; it is the 
usual coupling to gravitation which breaks the symmetry [4]. That model 
involves a principal bundle P ^ M with compact semisimple structure 
group 

G = SU(3) X SU(2) X U(l) 

as fiber; it postulates complex vector bundles associated to two representa- 
tions of G, whose sections are called Higgs and fermion fields. I'll leave the 

details of these representations unspecified (sec [5]): for our purposes the 
significant fact is that these are both bundles of (^n — l)-densities. 

The Lagrangian density of the standard model is a functional on a configu- 
ration space of fields (A, \I'), equal to the sum of 

• a Yang-Mills term defined by the curvature Fa of a connection 

one-form A on P, 

• a Dirac term ■ i^^^ = L{A, ^), and 

• a Higgs term L[A, $) of the form + -P(^)) with the latter term a 
polynomial in eg something like 

(with a dimensionless coupling constant A). 

The Yang-Mills term is conformally invariant in dimension four, and if we 
rewrite the auxiliary fields 

$ := ^ ,= *(n-2)/2n^^ 

in terms of fields ($o, ^o) of conformal weight zero, then the rescaling g 

L{A, *o) ^ L{A, *o), L{A, $o) ^ L{A, $o) 
with A:= A + u~^du. 

2.2.2 Sections g of the bundle 

G"^ XgP ^ M 

(where G'^'^ is defined by the conjugation action of G on itself) form a group 
9(G) of gauge transformations, which act on the space A of connections on 
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P by (g, ^) i-^ A — dg ■ g ^. The standard model Lagrangian is thus a 
function on the quotient 

(Higgs X Spinors) x^iG) 
(with g($, ^) = (g$, g^') on the Higgs and fermion fields). 

The analogous symmetry group for the geometric sector is the group D of 
(orientation and spin-structure-preserving) difFeomorphisms of M. This acts 
on S(G), and the moduli space of states for the standard model coupled to 
the usual version of general relativity is a bundle 

(Higgs X Spinors) Xg(G) ^1 ^ (• • • ) ^ rQ^/D . 

2.2.3 The multiplicative group rR(0)^ of nowhere- vanishing real- valued 
functions on M [§1.3] can be equally well regarded as a group S(Gm) of 
gauge transformations associated to a principal bundle with structure group 
the noncompact torus Gm(IR) = • If we interpret the fields of the standard 
model as densities as above, then its moduli space of field configurations can 
be presented as the quotient 

(Higgs X Spinors) Xg^g^ A 

with ^ 

A:=Ax dn\M), g(G) := g(G x Gm) = 9{G) x g((G^) 
under the action 

(g, u) ■ {A, *, $) := {A - u'^du - dg ■ g~^, u~^gW, u~^g^) . 

The Lagrangian of the standard model coupled to gravitation then extends, 
as in §2.1, to a function on the quotient of the space 

((Higgs X Spinors) x A) Xg(G^) (TQ; x r'R{^n - 1)) 

of fields by the gauge group S(G) x B. This defines a conformally invari- 
ant version of the standard model coupled to gravity, whose solutions are 
classical off the singular locus 7~^(0). 

2.2.4 This suggests several possible directions of extension, which I hope to 
discuss later: 

• the Connes-Chamseddine-Lott noncommutative version of the standard 
model [4] 

• supersymmetric (eg minimal) extensions of the standard model, 

• Connes-Kreimer-Marcolli renormalization of the conformally invariant 
4D model (which involves treating mass as a gauge field) [1]; and 

• questions of classical analysis: I've evaded certain issues by assuming 
the underlying space-time manifold to be compact, and by working entirely 
with smooth sections. However, §4.2 below suggests the interest of weak (eg 
Lf ) solutions to these equations. 
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§3 Inferences from scale invar iance 



" To show the fly the way out of the fly-bottle ..." 
Ludwig Wittengstein, Philosophical Investigations 

Invariance under rescaling by a constant factor p > may help illuminate 
some basic issues of physical interpretation: 

so 

2/n 



^ '-JM 



is scale-invariant, i.e. a 'pure number'. In view of the extensive literature 
concerned with anomalously large or small cosmological numbers, it seems 
remarkable that this L^-norm (when n = 4) seems to be well-behaved in 
the standard astrophysical models: the current bound for the cosmological 
constant is roughly 

|A| < 10-^^ sec-2 , 

while the universe is thought to be something like 4 x 10-*^^ seconds old, 
suggesting that 

\\R\\l2 ~ 0(1) 

(consistent with the hypothesis, plausible on other grounds, that R = 0). 

More generally, Holder's inequalities imply that for a smooth real-valued 
function / on M we have 

\\Ri9)f\\L^(g) < • ||/|liw(n-2)(,) . 

The norm appearing on the right rescales like a length: 



|^2n/(n-2)(-p4/(n-2)g) — /) 1 1 / 1 1 ^2n/ (n-2) (g) , 

which suggests regarding it as an estimate of the 'radius' of M (i.e. the nth 
root of its volume), measured in units defined by /. On the other hand, 
interpreting / as the inverse Planck length suggests regarding R^g)^ as an 
analog of the stress-energy scalar 

T ~ K-^R{g) . 

In fact 

ll^(5')/^llLi(p4/{n-2)g) = P^\\R{9)f\\L^{g) 

scales like ~ Energy^ • Hyper volume, cf. [8 §3.3]. The inequality 

\\R\\n/2 > P/'lll-|l/ll2"n/(n-2) 

can therefore be interpreted as bound of the form 

Mass 



Const > 



Radius ' 
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conceivably this lies behind the 'numerical coincidences' which physicists 
have interpreted as evidence for some version [3] of Mach's principle. 

§4 Two examples 

'"The further in you go, the bigger it gets," said Hannah 
Noon.' John Crowley, Little, Big 

I'll close with an attempt to show that this formalism is not completely 
without content. The first example below is quite widely known [cf. eg. [9 
§7]], but the second is more speculative. 

4.1.1 The positively curved Friedman pseudometric on = F_|_ is defined 
by 

' 1 
-R^A^l 

where A = (l + |r^)~"^ and = + + X3. The expansion factoi0 
R{t) {t = xq) satisfies the differential equation 

/dfi\2 Rq-R 



9 



— 



, i?(0) = 



\dtJ R 

it increases to a maximum Rq at time Iq, and then decreases to zero. The 
conformally equivalent pseudometric 

Rity^df - Si<fc<3 dxl 

on the stereographic completion M x S'^ is isometric to a product manifold 
with time parameter c?r = R{t)~^dt, i.e. such that 

{^f = {Ro-R)R, 

which is satisfied by R{t) = i?o sin^(r/2). In other words, the spherical 
Friedman model is conformally equivalent to the universal cover of the com- 
pact manifold x with the standard product Lorentz metric, endowed 
with a dilaton field (proportional to sin^(T/2)) which vanishes smoothly at 
singularities (big bangs) recurring with period 27r. [This corresponds to 
identifying the upper and lower edges in the Penrose diagram displayed in 
Fig. 21 (ii) of [8 §5.3].] 

4.1.2 The (hyper)volume of this model (i.e. of one cycle of the closed 
Friedman universe) seems not to be well known: it equals 

fto 



L 



*gl = 2-Yo\{S^)- r R{tfdt. 
Jo 



^I'm following the notation of [7 p. 117], but with ii denoting the Robertson- Walker 
factor. 
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If w := R/Rq then 



so 



and hence 



R^w' = (u;-^ - 1)-^/^ , 
dt = Ro{l - w)-^/'^w-^/^ dw 

r Ritfdt = Rl C nf'''^{l - w)-^''^ dw 
Jo Jo 



Substituting w = sin^ 9 yields 2 J^^"^ sin^ 9 dO for the right-hand integral; 
this equals 

— [60^ - 45 sin 20 + 9 sin W - sin 6^1 fJ'^ = — vr , 
96 ■'16 

yielding (tip o'the hat to Archimedes) the value 

2.2'K^.^'K.Rl = l-K^.Rt 
for the volume of one Friedman ^on. 

4.2.1 The Schwarzschild metric 

dr^ = (1 - 2mr-^)df - (1 - 2mr-'')-Ur'^ - r'^da^ 
is usually defined on the spacetime manifold R x — 0; 

dfj^ = + sin^ 0d02 
is the standard metric on the two-sphere. More generally, the expression 
dr^ = q{r)r~'^dt^ - q{r)~^r'^dr'^ - r'^da'^ 

with 

q{r) = — 2mr + = {r — m — D){r — m + D) 
defines the Reissner-Nordstrom metric; the discriminant = — will 
be assumed positive here. Both of these examples have vanishing scalar 
curvature, but if e > Reissner-Nordstrom space is not Ricci-flat. 

In null coordinates u, w such that t = ^{v + w) and 

\{v — w) = J q{r)~^r^dr 

the expression above defines the pseudometric 

dr^ = qr~^dvdw — r^da^ 

with associated volume element 

\g\'^dv Adw Ad9 Ad(j) = ^qfldv Adw = Ur^dr A dt , 

where Q = sin ^ d9 A d4> is the volume element for the two-sphere. The 
topology of the maximal analytic extension of such a pseudometric is com- 
plicated; see, for example, [8 §5.5]. Here we will be concerned mostly with 
the region of type II, in the terminology used there; this corresponds to the 
condition r G {m — D,m + D). 
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4.2.2 In these coordinates the Laplace-Beltrami operator 
takes the form 

A/ = 2g-n(rV,.),«, + (rV,^«),.] 
for a function / = f{v, w) independent of 9 and (j). If / is independent of t 
as weh, and we use primes to denote differentiation with respect to r, then 

v' = q~^r^ and 

A/ = 2q-\r^f'/q~'r'y/q-'r^ = 2T-\qf)' . 

There is thus a two-dimensional family of elementary harmonic functions 
u = u{r) in the interior Schwarzschild region, characterized by the condition 
u' = kq~^ for some constant k. Ignoring the time parameter, these functions 
arc all of Lebesgue class L^, but their derivatives are only locally of Lebesgue 
class L^. Since u" = —kq~'^q', any such solution has a point of inflection at 
r = m. The function defined by 

\ , , ,r — m — D . 

U{r) = l + -log — - 

m r — m + U 

if r G (m — D,m), and f/(r) = 1 if r > m, is particularly interesting. It is 
continuous, but not differentiable, at r = m; its derivative is 

U' = 2m'^D^q-^r] , 

where rj denotes a unit step function at r = m, so 

AU = 4m~^D^S 

with S a Dirac delta-function at r = m. The piecewise-differentiable tensor 
g = U^g thus has scalar curvature 

R = -24m-^D'^d 

which vanishes almost everywhere [14]. 

4.2.3 The conformally deformed tensor g will not be Ricci-fiat, even in the 
Schwarzschild case; a straightforward calculation [2 §6.3] shows that for a 
general conformal deformation g := 

Rl := RUg) = u-'/^-'RUg) + n-2"/"-2p^ 

with 

Pi = -2u[Vlu + (n - 2y^SiAu\ + 2(n - 2y^[nV'uVkU - 4|Vm|^] 

where signifies covariant differentiation [which on scalars is to be inter- 
preted as ordinary differentiation]. If we assume that n = 4 and u' = kq~^ as 
above, then the resulting tensor is diagonal in Schwarzschild coordinates; its 
entries have invariant significance, as the eigenvalues of considered as an 
endomorphism of the tangent space. Displaying these diagonal components 
as vectors, we find that equals 

kr-'^q-'^q'u[\, -1, 0, 0] + 2kr-^u[-\, -1,1,1]+ k'^r'^q'^ [1, -3, 1,1]; 
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if u' = kq^^rj is cut off at r = m, there is an additional term of the form 
—A;r~^(5[0, +1,0,0]. Assuming u = 1 at its inflection point, k = 2m~^D'^ is 
the unique value for which the determinant of vanishes at r = m; this 
characterizes the harmonic function U. 

It is similarly straightforward to show that (as in the Schwarzschild case), 
the equation of a radial timclikc geodesic in the metric g = U^g becomes 

U'^P = h^ + {2mr-^ - l)C/2 , 

where the dot denotes differentiation with respect to proper time, and h is 
a constant of integration. When r is small, 

is integrable; as in the classical case, such a geodesic reaches the origin in 
finite proper time. 

4.2.4 This suggests the interest of solutions of the equation 

Au + ku^ = 

with A 7^ 0. Taking e = and m = ^ for simplicity, and assuming as above 
that u depends only on the radial coordinate, this becomes 

2r-2(r(r - l)u')' + An^ = , 

which bears some (superficial?) resemblance to the Lane-Emden equation 
of astrophysics [6 Ch. IV]. 

By Cauchy-Kowalevskaya, the equation above has a unique solution analytic 
near of the form 

3(2A)-V2^(r) := 3{2A)-^/\-^/\l + Y,Wkr'') , 

k>0 

where 

(r(r - l)v'y + fr^-u^ = Q . 
In terms of v{r) = r~^/'^w{r), the equation above becomes 

4r^(r-l)u;" + 4r(2-r)u;' + 3(r-3)w; + 9w;^ = 0, 
which is satisfied by 

^26' W' 

Numerical computations (many thanks to S. Agarwala and Dan Christensen 
for invaluable help, including infinitely many corrections) suggest this series 
has radius of convergence one, and that w is nonvanishing in the interval 
[0,1). 

The (orientation-reversing) monotonic change of variables t = log|r~"^ — 1| 
maps (0, 1) to (—00,00). Rewriting the equation above as 
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suggests that v"{t) > for all t G M, and hence that the graph of v is concave 
upwards. Since v{r) — )• oo as r, s ^ 0, this implies the existence of a unique 
critical point v'{p) = with p e (0, 1). Numerical calculations suggest 

V{r) = r-3/2(l_^)-l/2(i_2^) 

as a reasonable approximation to v away from r = 1; it has a unique critical 
point at 

P = i(7- V13) - .85... . 
Reasoning as in §4.2.2, this suggests that the function 

V{r) = v{p)~^v{r) if r € (0,/?), = 1 otherwise, 

defines a conformal deformation of the interior Schwarzschild metric with a 
second-order phase transition at r = p, with 

R = QA = 108m^p-'^w{pf 

when r < p, such that 

^3/2 

w{p) 

has a logarithmic pole at r = 0, defining a complete metric which puts the 
singularity at infinity, possibly corresponding to an interesting new ground 
state for the interior Higgs field. The Penrose diagram for the associated 
conformally deformed metric glues together the vertical edges of Figure 25 
in [8 §5.5], identifying horizontal pairs of parallel type I regions. 

4.2.5 In terms of the coordinate t, the equation above has an asymptotic 
solution 

with coefficients in the differential Frechet algebra AC of smooth functions 
on M with rapidly decreasing (Schwartz class) derivatives. The function 

voit) = i[l + 8(1 + e-*)-2]V2 (= 1 [1 + 8(1 -r)2]V') G AC 

is an example: — >■ 1 as t — t- oo, as t — )■ — oo. If wc define 

v{r) := r-^/'^{l-r)-^/'^x{r) = cosh^{^t) x{t) 

then the equation for v can be rewritten as a Duffing equation 



c" + Six' + SqX 



2^ _ 9^3 



4 



X 



with 

_ 3e* + 2 - e-* ^ _ 9e* + 2 + e"* 
et + 2 + e-*''^0"4(e* + 2 + e-*) ^ 

This has xq = f^o (corresponding to fo) as a kind of asymptotically station- 
ary approximate solution. 
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To improve the approximation, let x = xoy (note that Sq is invertible in 
AC). The Hnear operator 

L := do'^id + 5i)d6o = dQ^^iSod^ + eid + eo) 

has coefficients 

ei = 2(5q + 5iSo, eo = Sq + SiSq 
in the Schwartz class S. It thus extends to define a map from the differential 
algebra AC[[t-i]] to §[[t-% 

Suppose now that y(n) = ^„>/j>o ?/fe*"'^ ^ AC[[i^^]] has been constructed, 
such that 

F(y(n)) := Ly(n) + y(n) - y(n)3 G r''-' AC[[t-']] ; 

we can start an induction with i/q = 1, since Ei = t{6Q + 6i6'q) G S. Then 

F(y(n)+y„+it-"-i) = F{y{n)) + F'{y{n)) ■ yr,+it~^~' mod t'^'^ AC[[t-']] 

with F{y{n)) = En+it~'^~^ modulo higher powers of t^^, for some asymp- 
totically constant error term E^+i- The coefficient of t"""^ in the term on 
the right above then simplifies to 

En+i + {L- 2)y„+i 

and we can take yn+i = (2 — L)~^-En+i- 
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